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PREFACE 


ro  S'  £>// 

This  report  is  a  generalization  and  extension  of  ORO-SP-24,  "  Optimum 
Dispersion  for  Gaussian  Salvo. "  The  first  paper  considered  circular  targets 
and  circularly  symmetrical  aiming  errors  and  salvo  dispersions  ( X  and  Y 
components  equal).  This  report  introduces  asymmetry  in  all  three:  elliptical 
targets  and  elliptical  Gaussian  distributions  of  aiming  error  and  salvo  dispersion. 

The  extension  permits  application  to  indirect  as  well  as  direct  fire  mis- 

i 

sions.  The  practical  applications  are  currently  under  study  and  will  be 
published  later,  together  with  a  nomogram  which  permits  ready  application 
of  the  method. 

Despite  the  general  interest  in  this  subject  and  the  many  closely  related 
papers,  the  authors  have  found  no  other  publication  which  provides  the  same 
results.  After  straightforward  mathematical  statement  of  the  problem,  this 
report  succeeds  by  approximations  in  collapsing  the  eight  major  variables  to 
a  smaller  number  (mostly  ratios)  which  are  amenable  to  practical  handling. 

This  report  is  separately  issued  to  permit  its  use  by  others  without 
delay,  and  to  generate  criticisms  which  might  be  useful  in  its  use  and  inter¬ 
pretation. 
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SUMMARY 


N  independent  sho^  each  of  kill  probability  p  and  ballistic  dispersions 

A 

cr  and  a  ,  are  fired  with  aiming  errors  t  and  t_,  at  an  elliptical  target  of 

A  y  x  y 

semi-axes  ax  and  ay.  The  probability  K,  that  the  target  is  killed  is  shown 
to  be  given  approximately  by 

r~. - —  ^  CO  -r  _  _ _ -u 


1  -  K  = 
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where  I  ftu)  is  a  Bessel  function, 


_  Np  axav 

“  ~  IT/ 

x  y 


1  /Sx»  V)  _  1  I  sx*  VI 

P-  2  +  fyli  ,1r~2l^-"7y2  I  ' 
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Equation  (19)  is  simplified  by  a  series  of  approximations  to: 


K  =  ©  C(Y,Z)  +  (1- 
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C(Y,  Z)  is  the  incomplete  gamma  function  and  <£(x)  is  the  central  area  of  the 

normal  distribution  (interpreted  as  zero  if  a  <  1).  In  the  instance  that  (j>  =  O, 


it  must  also  be  assumed  that  0  =  1.  0, 


♦Equation  numbers  refer  to  the  text 


CASUALTY  PROBABILITIES  OF  GAUSSIAN  SALVOS 

Consider  an  elliptical  target  with  semi  -axes  a  ^  and  ay .  Let 

the  aim-point  be  Gaussianly  distributed  about  the  target  center  with 

linear  standard  deviations  r  and  t  .  A  salvo  of  N  shot  is  fired,  each 

x  y 

of  which  has,  independent  of  the  rest,  a  conditional  probability  p  of  incapacitating 

the  target  if  it  hits.  Each  of  the  shot  is  Gaussianly  distributed  about  the  actual 

aim-point  with  standard  deviations  o  and  a  (so-called  ballistic  dispersions). 

x  y 

The  target  axes  and  the  principal  axes  of  the  two  Gaussians  are  all  assumed 

to  be  parallel.  The  problem  is  to  express  in  a  tractable  form  the  probability, 

K,  that  the  target  will  be  incapacitated.  A  simpler  version  of  this  problem 

was  treated  in  ORO-SP-24^/  with  the  restrictions  a  =  a  =  p,  <j  -  a  =  cr, 

x  y  x  y 

and  r  -  r  =  r.  The  probability  ,P(x,yj,  that  a  projectile  aimed  at  (x,y) 
x  y 

strikes  the  ellipse  may  be  written 
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To  obtain  an  approximation  to  the  above  integral  we  reviace  the  discrete 
target  by  a  diffuse  Gaussian  target.*^  We  write  the  Gaussian  target  as 


2tt  A  A 
x  y 
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(2) 
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and  equate  zero-  and  second-order  moments  with  those  of  the  ellipse.  Since 
the  Gaussian  target  distribution  is  centered  about  the  origin  the  zero-  and 
second-order  moments  are  clearly 


Poo  =  C,  P20  =  CAx2  and  p<,2  =  CAyl 


(3) 


For  the  ellipse,  let  x  =  a  r  cos  ©,  y  =  a  r  sin  ©,  dxdy  =  a  a  rdrd©  and 

x  y  x  y 

w’e  find 

Moo  =  7raxay »  Mm  *  rf-fady  =  aay  §  rdr  j 

ax  1  cos  W  0  ® 

O  oo 

=  -f-  ax»ay  and  simllarily  pj2  =  -|-axay3 

and  equating  the  two  sets  of  moments  (3)  and  (4)  yields 
C  =  7r a  a  j A  »  —  a  iA  *  -  a 

X  y  a  2  x.  y  jy 

and  (2)  becomes,  on  substituting  from  (5) 
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Thus,  for  the  sharp  target,  any  fragment  falling  within  the  ellipse  scores  1,  a 
fragment  falling  outside  scores  0y  while  for  the  diffuse  Gaussian  target  a  hit 
scores  2  at  the  origin  and  decreases  to  considerably  less  than  1  at  the  boundary 
of  the  ellipse  as  shown  below: 
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we  now  asK  for  the  approximation  to  the  probability  given  by  (1),  that  a  single 
shot  aimed  at  the  point  (x,y)  will  hit  the  ellipse.  Using  the  diffuse  target  (6), 


this  may  be  written 


P(x>y) 
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This  immediately  splits  into  the  product  of  two  integrals; 
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the  first  integral  can  be  written  as 
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Now,  it  is  well-known  that 
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md  the  first  integral  becomes 
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This  together  with  the  corresponding  expression  involving  y,  reduces  equation 
(7)  to 


P(x,y) 


{  a  a 
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Introducing  the  abbreviations  S  2  =  cr  2  +  4-  a  2  and  S  2 

xxx  y 

(8)  becomes  2 
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Having  assumed  independence,  the  conditional  probability,  K(x,  y),  that 
at  least  one  of  the  N  shot  incapacitates  the  target  is  given  by 

K(x,y)  =  1  -  (l-PP(x,yf  (10) 

with  p  the  conditional  probability  that  a  hit  will  be  a  casualty.  If  we  make 
use  of  the  Poisson  approximation 

1  —  (1  -  x)N  ss  1  -  e  ” ’^x  we  may  write  (10)  as 
-  Np  P(x,y) 

K(x,y)^l-e  .  (n) 

The  over-all  casualty  probability,  K,  is  obtained  by  averaging  (11)  over 
the  aim-point  distribution : 


Equations  (9),  (11)  and  (12)  now  yield 


1  — K; 
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where  we  have  written,  for  brevity  a  =  .  r  A-X 

2SxSy 


We  now  make  the  change  of  variable  x=  S  r  cos  ©  ,  y  —  S  r  sin  0, 

x  y 


dxdy  —  rdrdG  and  (13)  becomes 
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2ir  ~(H  00820  +^Sin2e)  2 
rdr  y  d  ©e  (14) 


One  thing  is  immediately  clear  from  (14);  namely,  that,  although  the  problem 

contains  eight  variables  (a  ,a  ,  cr  »<r  .  r  ,t  ,  N  and  p),  the  solution  involves 

xyxyxy  SS 

only  three,  which  might,  for  example,  be  taken  to  be  at  .  We 

x  Ty 

shall  soon  find,  however,  that  a  slightly  different  choice  for  the  latter  two 
variables  will  be  advantageous.  Putting  u  =  \  r2  in  (14)  we  find  that 
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Now,  for  the  Bessel  function  of  imaginary  argument,  we  hay©- 


1  ^  2t  1  z  cos  © 
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(16) 


Either  sign  in  the  exponent  may  be  taken;  in  other  words  Ig  is  an  even  function 
of  Z.  Therefore  we  may  write 


,  S  1  S  1 

and  we  see  at  once  that 
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Replacing  ©  by  0,  (15)  becomes 

«0 

1  -  K=f?  -  y  ■  f  duij 


4  x  ♦  y  u  cos  © 

e  d© .  (18) 


By  the  periodicity  of  the  cosine,  we  see  on  applying  (16)  to  (18)  that 
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This  is  the  desired  probability,  but  it  is  a  little  complicated. 

The  discussion  of  (19)  is  facilitated  by  introducing  an  auxiliary  quantity 
6  =  y/  .  From  (17)  it  is  clear  that 

6  5  7  *  P ,  hence  0  -  6  -  1  . 
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Indeed,  it  is  well  known-'  that 


z 

I<!  (z)  ~  e  (Z  —  «) 

/2~irz~~ 

and  consequently  (19)  does  not  converge  unless  y  /3.  The  case  6  =  0  has 
already  been  covered  by  a  contour  map.*^  One  possibility  for  handling  (19) 
might  be  to  have  it  computed  by  the  1103  and  similar  contour  maps  drawn 
for,  say  10  or  20  equally- spaced  values  of  6  between  0  and  1.  Such  a  port¬ 
folio  of  maps  would  not  be  particularly  difficult  to  work  with,  although  it 
would  in  general  be  necessary  to  read  on  the  maps  for  two  values  of  6 
(or  more, in  those  cases  where  K  is  a  sensitive  function  of  6)  and  interpo¬ 
late  to  the  exact  value.  The  chief  objection,  however, is  the  prohibitive 
amount  of  labor  required  to  draw  up  such  a  portfolio. 

And,  indeed,  an  over-scrupulous  insistence  on  a  precise  evaluation  of 
(19)  is  unwarranted,  since  (19)  itself  already  rests  on  the  diffuse  target  and 
Poisson  approximations  (cf.-i^ where  the  ranges  of  validity  of  these  are 
discussed).  To  approximate  (19)  by  a  simpler  expression  we  shall  first 
examine  the  behavior  at  6  -  0  and  fit  for  small  5  with  two  degrees  of  freedom. 
We  shall  find  this  fit  good  except  near  6*1.  We  shall  then  examine  the 
behavior  at  6  *  1  and  again  fit  with  two  degrees  of  freedom.  We  then  fit  for 
all  6  by  a  linear  combination  of  the  two  approximations,  thus  adding  another 
degree  of  freedom  to  our  fit.  The  result  is  an  expression  which  is  simple 
to  evaluate  and  which  is  highly  accurate  for  all  values  of  the  parameters. 

We  shall  require  certain  integrals: 
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The  roots  of  the  denominator  of  last  Integrand  are  given  by 


Now,  for  the  integrals  to  converge,  we  must  have  7  <  1,  so  that  both  ex¬ 
pressions  for  Z  are  real.  Since  their  product  equals  J,  one  lies  inside  the 
unit  circle  (is  less  than  1)  and  one  lies  outside  the  unit  circle  (greater  than 
1).  Only  the  former  contributes  its  residue  to  the  right-hand  side  of  (21). 
We  find  then. 
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Now  let  A.~..li  ~r.  =  a,  Z  =  W  +  A, 
y 


z  -  A-tY1  ~rl  =  w-  .AA  r.yl  =  w  -  b 
r  r 

and  the  residue  in  (22^  at  Z  -  A  becomes 
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and  we  obtain  finally 
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We  now  wish  to  Introduce  the  approximation 
*  fru)  De”du 


(24) 


To  determine  D  and  d  we  match  zero  and  first-order  moments.  Now  from  (23) 


with  k  »  0,  we  find 


)  du  Iq  (yu)  e~^a  . 


1  1 
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D®  du  s  T-  and  therefore  for  the  equality 


of  the  zeroth  moments 


By  similar  reasoning 


JW  -  ? 


( &  -  V) 


D 

=  <T"  * 


=  D 
^2- 


for  the  equality  of  the  first  moments.  Now  solving  (25)  and  (26)  for  d  and  D 
and  substituting  in  (24)  we  obtain 


h  (7*0  e 


-/Ju_  Sp- 
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Substituting  (27)  into  (19)  we  obtain  for  small  y 


l  -  k; 


^  du  e  -(^X^u  -  (y~, 


Let  us  investigate  how  well  (28)  may  be  expected  to  approximate  (19). 
When  r=  0  the  two  expressions  are  identical.  Also,  because  of  the  two 
degrees  of  freedom  in  the  fit,  one  anticipates  that,  as  y  increases  from  0  the 
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r* 


\ 

Thus,  both  solutions  approach  zero  as  y — >p,  but  quite  differently.  It  is  also 

easily  seen  that  the  values  of  K  given  by  (19)  and  (28)  both  approach  zero  as 

,  _ae-« 
a  ->  0  i  indeed  when  a  is  small,  the  factor  e  is  nearly  flat  for  all  u 

and  consequently  the  fit  we  have  made  in  (27),  which  is  perfect  for  a  constant 
factor  (a  =  0),  should  diverge  quite  slowly  frcm(26)as  a  increases  from  0. 

And  when  a  — both  expressions  for  K  approach  1. 

When  p  — >  0  the  fact  that  0  ~  at  -  p  shows  that  y  — >  0  also,  so  the  approxi¬ 
mation  is  automatically  good.  For  p  — >«> ,  by  replacing  u  by  ^  (19)  may  be 

written 

1-K.Xr?  ^  dulo(£  u)  e-'1-053 


Iq  is  an  even  function  cf  its  argument;  hence,  retaining  only  first-order  terms, 
we  get 


1  -  K 
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Likewise,  by  replacing  u  by  z  u  in  (28),  we  get 
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and  retaining  only  first-order  terms 
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Hence,  (19)  and  (28)  are  identical  to  the  first-order.  Thus  the  fit  is  also  very 
good  for  large  /3.  Therefore,  it  is  only  for  7-/?  that  we  need  to  expect  substan¬ 
tial  deviations  of  (28)  from  (19),  and  then  only  for  intermediate  values  of  a. 

1/ 

Reproduced  in  Figure  1  is  the  contour  map  given  already  in  ORO  SP-24 
for  the  function  we  will  denote  here  by 


c<y,z)  = 
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-Wu 


1 

2Y 


u 


(30) 


Comparison  of  equations  (28)  and  (30)  show  that  we  have  for  y^O 


K~C  •  2a  (//- 


(31) 


The  values  of  Y  and  Z  required  to  read  the  map  of  C  (  Y,  Z)  are,  recalling  that 
a  *  2|^Y  and  the  values  of  0  and  y  from  (17), 


Sx  Sy 

Npajjiay 


and  Z  - 


2Npaxay  S*' 


^x2 


(32) 


We  pass  on  now  to  the  fitting  of  (19)  in  the  case  where  y#p.  Before  we  can 
do  this  we  must  study  the  behavior  of  (19)  in  this  case.  We  have  already  seen 
that  K  — >  0  when  7  — >  (3  for  all  a.  However,  we  know  that  K  — >  1  as 
a  — > 00 ,  no  matter  what  are  the  values  of  /3  and  7.  Hence,  we  anticipate  a 
clash  when  a  is  large  and  7»  fl.  This  situation  is  pointed  up  in  (29),  where  the 
factor  in  front  of  the  sum  approaches  0  as  7 — >  (3,  while  the  sum  itself  (it 
is  not  difficult  to  show)  approaches as  a  — >».  One  might  be  tempted  to 
suggest  that  the  sum  in  (29)  be  plotted  up  as  a  contour  map  in  a  and  /3  and 
that  (29)  be  used  to  approximate  K  when  y~p.  This  would  mean,  however, 
that  by  taking  a  sufficiently  large,  K  would  be 
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estimated  to  be  bigger  than  unity,  and,  indeed,  as  large  as  one  pleases}  this 
is  not  the  sort  of  approximation  we  are  looking  for. 

What  is  happening  here  may  be  understood  by  visualizing  the  sequence  of 
contour  maps  as  7  — >  ft.  The  equiprobability  contours  keep  moving  out  into  the 
region  of  large  a.  Thus,  for  any  given  a.  these  contours  continue  to  slip 
past  as  7  — >  0,  so  that,  in  the  limit,  the  probability  equals  zero  for  all  finite 
a.  In  order  to  stay  in  the  region  of  non-zero  probabilities,  we  must  let  a 
become  large  at  the  same  time  that  we  let  7  — >  ft.  This  is  something  we 
would  do  automatically  if  we  were  drawing  up  a  portfolio  of  contour  maps, 
since  we  would  be  guided  not  by  the  magnitude  of  a  but  by  the  desire  to  in¬ 
clude  all  contours  of  reasonable  probabilities. 

Now  the  function  f(u)  5  e  *"  ae  ,  which  occurs  in  (19)  has  a  maximum ; 
e~p(^)^  for  u  s  log  ~  .  If  we  employ  the  fact  that  Iq  (z  — >°)  in 

(19)  we  find  the  integrand  of  (19)  • 


e~(ft~y) u  ae  *  jf  we  ignore 


the  slowly- varying  factor,  -fir ,  here,  which  cannot  shift  the  maximum  by  much, 

we  find  that  the  maximum  of  the  integrand  occurs  at  u  z  log-^ — ^  Hence 

when  yzfi,  this  maximum  will  occur  at  a  large  value  of  u  and  thus  only  large 

u  will  contribute  appreciably  to  the  integral  in  (19).  Therefore,  the  closer 

z 


Y  is  to  0  the  better  the  asymptdtic  expression  - 


e 


approximates  the 


Bessel  function  so  far  as  the  integration  in  (19)  is  concerned.  Making  this 


approximation  in  (19)  we  have  as  7 


(33) 
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(We  see  incidentally  from  this  why  (28)  does  not  reproduce  the  behavior  of 
(19)  accurately  when  y^p.  For  the  behavior  of  e-^,  (yu)  is  dominated  by  the 
exponential  factor  e  '  u  for  large  u,  so  long  as  y  differs  appreciably  from 
Pi  and  hence  (24)  is  a  good  approximation  to  useu  When,  however,  yzP,  the 
exponential  factor  becomes  of  correspondingly  less  importance;  the  integra¬ 
tion  in  (19)  begins  to  be  sensitive  to  the  square  root  factor;  and  (24)  is  no 
longer  a  good  approximation  .) 

Writing  u  =  Z2  in  (33)  we  find 

1  -  \  e-(0-7)  Z  -ae  2dZ,  Qr 

o 

1-K^=5_  J  zJ-°»  dZ  t(y—*P),  (34) 

— «o 

since  the  integrand  is  now  an  even  function  of  Z*  In  other  words, 

-oe“z2  1 

1  -  K  ~<  e  >  where  Z  =  N  (0,  )»  (35) 


where  the  corners  indicate  the  expected  value, and  the  latter  half  of  the  equation 

states  that  Z  is  normally  distributed  with  zero  mean  and  variance  ^ 

From  this  we  see  that  when  the  variance  of  Z  is  large  and  so  the  ex- 

pection  in  (35)  gives  more  and  more  weight  to  the  wings  of  the  function 
-Z2 

e  ,  i.e.,  since  the  latter  approaches  1  as  Z  — this  says  that,  for 
finite  a,  1  -  K  — >  1  as  7  — >  p,  verifying  what  we  noticed  before. 
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It  is  convenient  to  rewrite  (35)  as 
K~<  I-.-8"*  e-^>Z!  dZ  .  (36, 

If  in  the  integral,  we  set?  =  0,  we  find  that  it  still  converges,  call  it 
~  {  „-Z*\ 

J(a)£  r  •  ■ 1 

At  Z  s  0,  the  integrand  equals  1  -  e“a,  and  a  'typical'  value  of  the  integrand, 
somewhat  analogous  to  the  half  life  of  a  radioactive  substance, will  be  gotten 
by  setting  Z  =VI5il T.  Thus,  for  large  a,  we  may  expect  the  area  under  the 
curve  to  be  of  the  order  of  the  magnitude  of  the  height,  1  •  e  times  the 
typical  width,  Vlogar  x 

J  (a)/v/A^log  a  (a — **)  • 

Examination  of  the  previous  three  equations  shows  that  real  and  non-zero 
values  of  K  are  obtained  only  when 

a  s  V3-7  (o  >  1  and  rj  >  1)  ,  (38> 

where  tj  is  some  new  finite  value.  This  is  precisely  the  sort  of  behavior 
anticipated,  for,  when  rj>  1,  a  «  as  y — >pt  which  corresponds  to  the 
upward  shifting  of  the  probability  contours  described  earlier. 

Making  the  change  of  variable 


y? 
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in  (36),  and  using  (38),  we  obtain 

~  4f£ 

L  (1 '  "Y 


K'v ' 


)  e-yS  l0®  **  dy  . 


Now,  when  tj  >1,  as  we  are  assuming  here. 


1-e"* 


1~ya 

P~7 


-*  i  (y*<  i) 

->  o  (y*  >  i) 


as  y  • 


B 


Thus 


-y2  logT? 


dy^(7 — >0)  . 


(39) 


Defining  the  cumulative  normal  distribution  as  usual  by 


fe-  c  « 


?  w  *  vs? 

(39)  may  be  written 

( v/2  log  tj  )  (tj>l) 


as  7— */3  . 


—  °  (tj<  1)  , 

Finally,  recalling  that  tj  =  our  result  becomes 

K~  ^  (V/2  (j8-t)  log  a  ^  (or  >1) 

~-'0  (cr<I) 


as  7 — >/3  . 


(40) 


(41) 


(42) 


It  is  now  a  very  simple  matter  to  improve  our  result  by  fitting  (19) 
to  it,  analogously  to  the  way  we  have  already  handled  the  case  7~0.  We  want 


-flu  „  JL  -hu 

e  4  (>u)  <•  ® 


(43) 


18 


Solving  for  H  and  h,  we  get 


Substituting  in  (43)  gives 


and  substituting  (44)  in  (19),  we  find 


and  again  writing  u  =  Z2 


Hence  we  can  write  the  solution  from  (42) 


Finally,  it  remains  to  splice  togther  the  two  solutions  (28)  and  (46). 
To  do  so,  we  now  make  the  approximation 

(*yu)  **  a"  U  +  (!-&)  — —  e"  u  (47) 

P  V2irj8u 

i.e.,  we  make  a  linear  combination  of  equations- (27)  and  (44).  Equality 
of  the  zero  and  first-order  moments  of  both  sides  is  guaranteed  by  the 
form  of  the  expression  on  the  right.  Thus,  to  determine  ©,  we  match  the 
second  moments  of  both  sides.  We  now  use  (23)  with  k  =  2  and  find 

W'+v'.-e  hi  -  9|  — M!_ ^ 

(13-  -  y2>  2  (ft’  -  y2)*5  (ft’-  -  y2)'5  • 

from  which  it  follows  immediately  that 


Substituting  (47)  and  (48)  into  (19)  and  writing  the  answer  in  terms  of  (31) 
and  (46)  we  obtain  our  final  result 

K  «(1  -y)C  •2a(^-V))  ^  ^  )  lo&  a*  (49> 

where  ^  is  to  be  interpreted  as  zero  if  a<l.  We  now  recall  the  definitions 

of  quantities  occuring  in  the  equation  in  terms  of  the  parameters  of  the 
original  formulation  as  follows: 


a 


ay> 


P 


,  S  2  S  2 
1-  x  y 

7  =  2"  v"  "  V 


(50) 
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| 


JLt 


should  b©  recalled 


Equation  (32),  which  tells  us  where  to  read  on  Figure  1  is  still  valid;  and 

-j  v 

we  note  also  that  the  quantity  0  =  1-  (-s-)2  which  occurs  in  two  places  in 

P 

(49>l,mere.ye  =  l-£{|(^+^)]  .It 

that  <j>  Is  to  be  interpreted  as  zero  if  a  as  1.  In  this  instance  0  must  also  be 
interpreted  as  equal  to  1.  0.  ©  is  defined  as  a  weighting  factor  between  two 

solutions,  if  one  solution  does  not  apply,  the  full  value  of  the  other  solution 
must  be  used. 

From  the  form  of  the  weighting  coefficients  in  (49)  we  see  that  for  ^ 
small,  the  second  term  contributes  very  little  to  the  answer;  indeed  the 
coefficients  have  a  horizontal  tangent  at  y  =  0.  In  the  latter  case,  however, 
when  a  slight  drop  in  ^  produces  a  proportional  admixture  of  the 

first  term  to  the  second  term.  In  other  words  (31)  remains  accurate  over 
a  wider  region  than  (46). 

From  the  form  of  (38)  we  see  that  when  ysP  and  a  approaches  1, 
then  tj  approaches  1,  but  still  more  rapidly.  Since  the  validity  of  the 
function  <j)  depends  on  tj  being  more  than  I,  it  is  to  be  expected  that  the 
solution  has  a  danger  point  that  should  not  be  closely  approached . 

This  turns  out  to  be  the  case,  and  the  critical  region  appears  when  <|) 
exists  and  a  is  approximately  1.  This  situation  is  easily  handled  in 
practice  by  raising  a  somewhat,  perhaps  by  increasing  N,  and  then 
extrapolating  back  to  the  desired  value  of  a. 
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An  illustrative  example  of  the  application  of  the  formulas  developed 
in  this  paper  is  included  in  these  final  paragraphs.  The  particular  formulas 
used  in  any  application  are  scattered  throughout  the  paper.  The  relations 

between  S2 ,  o2,  and  a2  were  introduced  immediately  following  equation  (8). 

a,  as  a  function  of  N,  p,  a  ,  a  ,  S  «  and  S  ,  was  first  introduced  following 

x  y  x  y 

equation  (13)  and  restated  in  equation  (50).  /3  and  y  are  shown  in  equations 

(17)  and  (50).  The  C(Y,  Z)  function  is  graphed  in  Figure  1;  both  Y  and  Z  are  defined 
in  equation  (32).  One  now  requires  only  0  and  <£)  to  determine  the  probability  of 
a  casualty,  K.  0  is  defined  in  equation  (48),  and  <£>  is  the  central  area  of  the 
Normal  (Gaussian)  distribution.  The  final  probability  of  a  casualty,  K,  is 
determined  in  equation  (49). 

A  Ritnnle  example  is  selected  which  affords  a  comparison  with  earlier 

work:  A  weapon  that  has  a  large  dispersion  characterized  by  range  and  deflection 

probable  errors  of  94  and  18  yards  respectively  is  used  to  attack  a  small  target 

4  yards  in  radius  (symmetrical  target,  a  =  a  ).  A  single  hit  destroys  the 

x  y 

target  (p  =  1).  What  is  the  probability  of  destruction  for  32  shots  for  aiming 

error  CEP's  of  0,  50,  100,  and  200  yards  (symmetrical  aiming  error,  r  =  r  )  ? 

x  y 

In  converting  from  range  and  deflection  probable  errors  to  standard  deviations, 

one  obtains  tr  =  140  yards  and  a  -  27  yards.  As  a  result 
x  y 

S2  =  19,  600 
x  ’ 

and 

S2  =  733. 

y 
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By  the  relations  given  in  equation  (17)  and  (50) 

0=  0.0677 


y  =  — —  =7.38. 

20! 

The  relations  involving  ft  and  y  both  depend  on  the  evaluation  of  r,  given  in 

this  example  In  terms  of  CEP's.  Because  of  the  symmetry,  r=  r  =  r  . 

_ x  y 

Once  r  is  determined,  ft,  y,  Z,  O,  and  | P ©log^  a  can  be  obtained.  The 
following  table  shows  these  figures  : 

Table  1 

PARAMETERS  VS  CEP 


CEP 

T 

P 

y 

z 

e 

^©logea 

0 

0 

1 

- 

- 

- 

- 

50 

42 

5.64 

5.23 

4.74 

0.139 

-2.05 

100 

85 

1.41 

1.31 

18.8 

0. 139 

-0.511 

200 

170 

0.352 

0.327 

75.8 

0.139 

-0. 128 

300 

255 

0. 157 

0.145 

171 

0.139 

-0. 057 

In  each  case  vJ/ieiog~a  is  imaginary,  and  we  must  take  <j>  =  0.  When 
,J”j9  ©loge<*  is  positive,  ^  is  just  the  central  area  of  Normal  (Gaussian) 
distribution.  When  <[>  is  zero,  0  must  be  assumed  equal  to  1. 00.  In  this 
example  (j)  is  always  zero,  0  is  always  1. 00,  and  X,  the  probability  of  a 
casualty,  is  just  C(Y,  Z).  The  results  are  given  in  the  following  table: 


T 


